Anisotropy of graphite optical conductivity 
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The graphite conductivity is evaluated for frequencies between 0.1 eV, the energy of the order of 
the electron-hole overlap, and 1.5 eV, the electron nearest hopping energy. The in-plane conductivity 
per single atomic sheet is close to the universal graphene conductivity e^ /Ah and, however, contains 
a singularity conditioned by peculiarities of the electron dispersion. The conductivity is less in the 
c— direction by the factor of the order of 0.01 governed by electron hopping in this direction. 

PACS numbers: 78.67.-n, 81.05.Bx, 81.05. Uw 
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Recently, the light transmittance of graphene was 
found [iHsl in the wide frequency region to differ from 
unity by the value of ira, where a is the fine structure 
constant of quantum electrodynamics. These experimen- 
tal observations are in excellent agreement with the the- 
oretical calculations d, Q of the graphene conductance, 
G = /Afi, which does not depend on any material pa- 
rameters. 

This phenomenon is remarkable in two aspects. First, 
the fine structure constant has been found in one mea- 
surement for the first time in solid state physics. Second 
and most important, the Coulomb interaction does not 
disturb the agreement between the experiment and the 
theory d, 0| • It should be emphasize that the Coulomb 
interaction in graphene is poorly screened while the car- 
riers are absent in this gapless insulator. 

In connection with this, it is interesting to study 
the change in the optical conductivity going from 2d 
graphene to its close "relative" , 3d graphite, with the 
optical conductivity measured in Refs. 0,0. 

The electron properties of graphite is well described 
within the classical Walles-Slonczcwski-Weiss-McClure 
theory [loj . There are many parameters in this theory of 
the various order of value (see, e. g. [lH). Among them, 
the energy 70 = 3.1 eV is largest one representing the 
electron in-plane hopping between nearest neighbors at 
the distance rp =1.42 A. If we are interested in frequen- 
cies less than 3.1 cV, we can use the power k— momentum 
expansion of the corresponding term in the Hamiltonian, 
taking only the linear approximation. Then the constant 
velocity v = 10® cm/s appears. The parameter 71^ 0.4 
eV known from optical studies of bilayer graphene jl^lTsj 
is next in the order. It describes the interaction between 
the nearest layers at the distance co=3.35 A. The pa- 
rameters 73 and 74 give corrections of the order of 10% 
to the velocity v. Finally, the parameters 72, 75 of the 
order of 0.02 cV from the third sphere arc used in order 
to describe the dispersion of the conduction and valence 
bands in the c— direction. They arc usually included in 
order to characterize the carriers and are known from the 
Shubnikov-de Haas oscillations and the cyclotron reso- 
nance. However, for the optical transitions at relative 
high frequencies 72, 75 <C a; <C 70, we can, first, neglect 




FIG. 1: Dispersion of low energy bands in graphite. 



the smallest parameters 72,75 and, second, use the lin- 
ear k— expansion with the constant velocity v for in-layer 
directions. Our results have the explicit analytic form. 
Thus, the simplified Hamiltonian of the model is given 

by 
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where the velocity parameter v is included in the defini- 
tion of the momentum components k± — v{^ikx — ky), 
and the constant 71 stands in the function 7(2) ~ 
271 cosz depending on the dimensionless fc^— component 
z = k^c^ along the c-axis, < z < 7r/2. 

The eigenenergies of the Hamiltonian write: 



Si2 = ^±J^+fc^. 



(2) 



e34 = -4^±J^ + fc^. 
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The so-called "Dirac" point of graphene, k = 0, turns 
into the K-G-H line of the graphite Brillouin zone, where 
the valence and conduction bands slick together, £2,3 — 0. 
It should be emphasized that this degeneration is condi- 
tioned by the lattice symmetry but is not a result of the 
model. 

Others two bands, £14 = ±7(2), are spaced at the 
distance 7(2) which vanishes at the H point of the Bril- 
louin zone. This band schema corresponds to the gapless 
semiconductor. 

In order to calculate the optical conductivity, we use 
the general expression Q 
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valid in the collisionless limit lo ^ v, where v is the col- 
lision rate. This condition is definitely fulfilled, if the 
frequencies are larger than the electron-hole overlap in 
graphite determined by the parameters 72 , 75 • The tem- 
perature is involved here by the Fermi-Dirac function 
/(er) = [exp(^^) -I- 1]""^, the coefficient 2 takes into ac- 
count the spin degeneration, and the integral is taken 
over the Brillouin zone where the electron dispersions £„ 
are defined. 

The first term in Eq. ([3]) is the intraband Drudc- 
Boltzmann conductivity with the group velocity 

v„„ = dsn/d'k. 

This conductivity behaves as 1 /w and becomes less than 
the second term for frequencies higher than the electron- 
hole overlap. The second term corresponds with the elec- 
tronic interband transitions accompanied by the photon 
absorption. It involves the matrix elements of the veloc- 
ity operator 

OK 

calculated in the representation transforming the Hamil- 
tonian ^ to the diagonal form with the help of the op- 
erator U . We find for various transitions 

2i(£3 - £2)ky/N2N3 , 
--2{ei+e2)kjNiN2, 
2i(£4 - ei)ky/NiN4 , 

where = 2{el + fc^) . 

The calculations show that the off-diagonal compo- 
nents of conductivity reduce to zero and the in-plane 
diagonal components are equal. For their real part, we 
obtain the integral which is explicitly taken over ip and k 
in the polar coordinates at the zero temperature. Thus, 
we meet the integral over k^ : 

<.7r/2 
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where 7(2:) = 271 cosz, and 0i, O2 are the step functions 
depending on w — 7(2;) and w — 27(2), respectively. This 
integral can also be taken, but the result looks more com- 
plicated. 

Here, we introduce the conductivity ctq = e^/47ico 
which can be named as the graphite universal conduc- 
tivity. It differs from the graphene conductivity only in 
the factor I/cq which is simply the number of the atomic 
sheets in graphite per length unit in the c— direction. One 
can see, that the graphite conductivity goes to ctq at low 
as well as high frequencies compared to 7(0) (see, also 
Fig. (2). However, at w = 271=0. 84 eV, both the kink and 
the threshold are seen in the real and imaginary parts, 
correspondingly. These singularities arise due to the elec- 
tron transitions between bands 2^1 and 4 — 3 (see. 
Fig. [T]) described by the second term in Eq. (j?]). The 
position of the singularities gives the value 71=0. 42 eV, 
which agrees well with optical studies of bilayer graphene. 
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FIG. 2: Real cti and imaginary 02 parts of the graphite optical 
conductivity for the in-plane direction (per one atomic sheet 
in units of /Ah) versus frequency (in units of 271 



eV); experimental data 
theory, dashed lines. 
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Let us consider next the conductivity in the c— axis. 
We need now the matrix elements v^^. Calculations 
show that they are nonzero only for the transitions 2 — 1 
and 4 -)> 3: 

v'21 = 27'(z)£i£2/7Vi7V2 , 

Vl^ = -27'(z)£3£4/Af3A^4 , 

where the derivative ^'(z) = 27icosin2:. Using Eq. ([2]), 
we get 



-vl^ = -7'(2)fc/v/72(z)+4fc2. 
Integrating in Eq. (jS]) over Lp and fc, we obtain 
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FIG. 3: The real and imaginary parts of conductivity in 
c— direction; units are the same as in Fig. [2l 

where the integral over 

4 r^"^ ( cos^z\ 
I{t) = -J fizsiii^z|^l — je{t- cos z) 

with t = a;/27i. This integral can be also taken and it 
has in limiting cases the very simple forms: 
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The imaginary part of the conductivity in c— direction 
is given by the fc^— integral 
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The conductivity in the direction is shown in 
Fig. 131 Compared with the in-plane conductivity, the 
c— conductivity is less by the factor (7iCo/?iw)^ ~ 0.01. 
This factor represents the squared ratio of the hopping 
integrals for the inter- and in-layer directions (71/70)^ — 
exp(-2co/ro). 

In conclusions, for the in-plane direction, the optical 
conductivity of graphite per single atomic sheet is close 
to the graphene universal conductivity. However, the sin- 
gularities, the kink in the real part and the threshold in 
the imaginary part, appear at the frequency uj = 271, 
where 71 is the inter-layer hopping energy for the bilayer 
graphene. For the c— direction, the conductivity is less 
by the parameter representing the ratio of the inter- and 
in-layer hopping energies; the real part of conductivity 
increases linearly with the frequency and does not con- 
tain any singularities. 
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